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Abstract 

In this paper we study the spectra of glueballs on the Klebanov-Strassler back- 
ground and its extension to the baryonic branch. We numericahy calculate the mass 
spectrum of glueballs from the spin 2 "gravity" multiplet, which contains the trace- 
less part of the stress-energy tensor and the transverse part of the U{1) 7^-current. 
The mass spectra of the corresponding fluctuations in supergravity coincide due to 
supersymmetry, which is manifest in the effective five-dimensional theory through 
a Supersymmetric Quantum Mechanics transformation. We show that the glueball 
spectra grow as ccUn? for large values of the baryonic branch parameter U . 



1 Introduction 



Ideas of holography underlaid the AdS/CFT correspondence [T] provide a promising per- 
spective for the study of gauge theories via string theory and supergravity. In particular, 
the extension of the holographic principle to a non-conformal case enables to capture the 
strongly coupled dynamics of a gauge theory through the classical supergravity. This al- 
lows, at least in principle, to calculate various correlation functions, extracting the masses 
of glueballs, which is not possible by means of the standard field theory technique. In 
this paper we focus on a particular case of the Klebanov-Strassler supergravity back- 
ground [2] and its extension to the baryonic branch [21 SJ [5] , which is a gravity dual of 
the non-conformal Af = 1 gauge theory with two pairs of matter multiplets. 

According to the AdS/CFT correspondence, gauge theory operators correspond to 
fluctuations of the background supergravity fields. Thus the stress-energy tensor cor- 
responds to the fluctuation of the metric. The transverse traceless part of the former 
combines with the transverse part of the U{l)n current and the transverse fermionic 
superconformal current into a spin 2 massive super mult iplet [6] . The mass spectra of the 
corresponding glueballs coincide, what is evident from the supersymmetric structure of 
the equations of motion in the gravity dual theory. 

In the gauge theory the supersymmetric structure of the spin 2 multiplet is transparent 
through the on-shell equation of current conservation [6j . Using superfield notations, one 
has 

^"Ka = D^S, (1) 
where Vaa = ^^acXtJ^ is ^ ^^^^ supercurrent that contains the T'^" and the J5 current, 

= J^-l^ Q^y"^ T/ + ie^d^'s-i P + i Pe"^ (2D^ + UJ^) + fermions. (2) 

The supercurrent contains two supermultiplets: the transverse spin-2 multiplet, which 
consists of the traceless transverse components of T^^ and the transverse part of J^, and 
the chiral multiplet S containing the trace T^^, the divergence 9^J^ and the 7-trace 
of the superconformal current. It also contains the complex scalar field s as its lowest 
component. The chiral multiplet 5* thus accounts for the anomalies of the scale, U{1)ti 
and the superconformal symmetries, associated with the components of the supercurrent, 
while the equation ([T]) is the supersymmetric generalization of the anomalous divergence 
of the current. 

On the gravity dual side, background fluctuations of supergravity fields that are dual to 
the operators in ([2]) are massless in the five-dimensional sense if the theory is conformal. If 
the conformal symmetry is broken, the dual of the U{l)'ji current satisfies the equation for 
a massive vector particle in five dimensions. Note that the dual of the traceless transverse 
part of the stress-energy tensor is described by a five- dimensional massless equations in 
both cases. However the finite warp-factor at the tip of the conifold in the non-conformal 
case leads to a finite four-dimensional spectrum of glueballs. 

In this work we describe the holographic dual modes of the traceless part of ([2]), 
namely the spin 2 gravity multiplet, in the context of the baryonic branch of the KS 
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background. The baryonic branch is a continuous family of the type IIB supergravity 
solutions originating at the KS background. The branch is parameterized by the vevs of 
the baryonic operators in the dual gauge theory [HITj. The backgrounds from the family 
are constructed in terms of the Papadopoulos and Tseytlin [S] ansatz, which consists of 
scalar functions parameterizing the metric and fluxes. Those scalar functions depend only 
on the radial coordinate of the conifold t and satisfy a system of first order differential 
equations, which was derived in [1]. No analytical solution to this system is known, except 
for the KS cas^, and in practice the backgrounds from the baryonic branch have to be 
constructed numerically. More details about the baryonic branch solutions can be found 
in the works [U [5] . 

The KS and baryonic branch backgrounds correspond to non-conformal gauge theories. 
We derive the linearized equations for the vector fluctuation dual to the U{1) 7^-current 
and the metric fluctuation dual to the stress-energy tensor. Then we numerically compute 
the four-dimensional mass spectra of the corresponding glueballs along the branch. Since 
we derive our equations on the solid ground of "microscopic" ten-dimensional theory, 
we can test the applicability of the similar results obtained in [10] through an effective 
approach of five-dimensional models of gauge/gravity correspondence. 

This paper is organized as follows. In the section [2] we remind the reader the dual 
description of the bosonic operators of the gravity and anomaly multiplets in the case 
of the KS background. This part also contains a sketch of the holographic anomaly 
mechanism suggested in [TT] . 

Sections [3] and m are dedicated to a derivation of the linearized equations for the bosonic 
fluctuations of fields dual to the gravity multiplet of ([2]). On the supergravity side the 
bosonic sector of the multiplet consists of the symmetric traceless perturbation of the 
metric - the graviton, and the transverse vector perturbation discussed in the section |5J 

Since the bosonic fluctuations of the gravity multiplet are related by supersymmetry, 
their spectra coincide. There is a Supersymmetric Quantum Mechanics (SQM) transfor- 
mation relating the effective five- dimensional equations, which is a reminiscence of the 
original supergravity transformation in ten dimensions. In section H] we derive the equa- 
tion for the vector mode only for the case of the KS background, but the supersymmetric 
structure of the equations allows us to extend it further to the baryonic branch. We show 
that this equation is the same as discovered by [TO] in the five-dimensional approach. This 
is discussed in detail in section [51 

We present the results of a numerical calculation of the spectrum in the section [61 
Although the equations that describe the graviton and the vector particle yield the same 
spectrum of bound states, they are essentially different. We perform two separate calcu- 
lations of the spectrum of the gravity multiplet which is an important consistency check 
of the numerical results. We conclude with a discussion in section [71 

-'^The analytical solution f9], known as Maldacena-Nunez background, also solves the system of [4]. 
However, it has different boundary conditions at infinity and therefore does not belong to the baryonic 
branch. 
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2 Multiplets and Anomalies in the Dual Theory 



The purpose of this work is the study of the gravity multiplet, i.e. the fluctuations above 
a classical supergravity background dual to a field theory supermultiplet consisting of the 
traceless part of the stress-energy tensor T'^'^, spin 3/2 conformal supercurrent and the 
conserved part of the U{l)-ji current J5 . The gravity multiplet therefore contains the 
traceless symmetric excitation of the metric - the graviton hfj_^, the spin 3/2 gravitino 
and the transverse vector excitation along the 1-form dip which we specify below. 
Classical supergravity backgrounds that we study in this work are the backgrounds from 
the baryonic branch of the KS solution. 

First we consider the transverse non-diagonal fluctuation of the background metric 

9^iu — ~\~ h^i/, (3) 

where has only components in the Minkowskian directions. We find that in accordance 
with the general results of [12], [13] this excitation is described by the massless scalar 
minimally coupled to the metric in the Einstein frame for all backgrounds of the baryonic 
branch. 

Next we consider fluctuation that is dual to the conserved (transverse) part of the 
U{l)'ji current in the KS background. Recall that the KS solution corresponds to a 
manifold that is locally a product of the Minkowski space-time and the six dimensional 
deformed conifold [2], 

dsls = h-^/\t)r]^,dx^dx'' + h'/\t)dsl (4) 

where h{t) is the warp factor that depends on the radial coordinate t of the conifold, related 
to the standard conical radial coordinate r via t ~ 31ogr. U{1)ti transformations act 
as rotations along the conifold base T^'^, 

where %jj is one of the angles on the baseH 

In the conformal case, the background is invariant under this symmetry, what results 
in a massless gauge field A^. The KS background, as well as the backgrounds along the 
baryonic branch, breaks the U{1)ti symmetry already in the UV. The 2- form potential 
for the RR form F3 has an explicit ip dependence. In the UV limit 

where M is the flux of F3 through the of T^'^, and uj2 is the ip independent 2-form on 
T^'^. Given that ip itself is a double cover of the circle, C2 breaks U{1)ti down to Z2M in 
the UV. In the IR the metric has an explicit iJj dependence that breaks Z2M further to Z2 
in the full agreement with the gauge theory. 

^One can think of T^'^ as of the space x S^/U{1) ~ x S^. The angle ip is obtained by the 
identification of the 3rd Euler angles of both . 
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As a result, the corresponding fluctuation of the background acquires mass that is not 
vanishing even in the UV region [Ullllj. The fluctuation in question modifies the metric 
along the ip direction (7^^ and can be described by the perturbation of the 1-form dip by 
the "gauge" field A = A,, dx^ + At dt, 

dip ^ dip + A. (5) 

Since the dependence on the angles of the conifold is not important, we can restrict 
our attention to the five- dimensional theory. In the conformal case, in the absence of the 
3-form fluxes, the five-dimensional vector field A satisfies the equation for the massless 
vector 

d*5di = 0. (6) 

The longitudinal part of A is not fixed by the equation as it is a gauge degree of 
freedom. The corresponding symmetry is anomaly free. After adding the fluxes, the 
equation for A can be brought to the form 

d(/*5dM)) + *5i = 0, (7) 

with some background-dependent functions / and g. The longitudinal part of A is no 
longer trivial and satisfies 

d*5i = 0. (8) 

For an observer in four dimensions, the five-dimensional no-source equation is precisely 
the equation with an anomalous source 

d^A, = e{A) , (9) 

where n denotes the space-time indices. This holographic anomaly mechanism is discussed 
in more detail in [llj. 

The backgrounds we are interested in have a global SU (2) x SU (2) symmetry. Since 
we are interested in the uncharged sector, all fluctuations should be s-waves with respect 
to the directions along the base of the conifold. This is obvious for the four- dimensional 
metric fluctuations as we keep it angle- independent. In the case of the vector, it is more 
tricky. In fact we need to switch from the 1-form dip to the invariant extension —>■ g^+A, 
where 

= dlp + COS 6'id0i + COS 6'2d02 . (10) 

Apparently the shift of dip results in the same shift of g^. 

The anomalies of the scale, superconformal and U{1)ti symmetries form a chiral super- 
multiplet [B]. Its bosonic part on the gravity side contains the fluctuations of the metric 
trace hj^" and the longitudinal part of the vector field = d^a. 

In the section H] we derive the equation ([7]) for the transverse part of the vector fluc- 
tuation ([5]). The transverse component decouples from the longitudinal part and from 
other supergravity fluctuations. Unfortunately it is much more complicated to derive the 
equation for the longitudinal mode d. Moreover there are certain indications that it does 
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not decouple from the other fields and needs to be considered as a part of a more com- 
plicated system p5]. Coupling with different supergravity excitations will lead to some 
non-trivial right hand side of the equation iQ. It is particularly interesting to find the 
supergravity expression for 6 (A) and compare it with the gauge theory predictions. This 
task is more ambitious and we leave it for a future work. 



3 Graviton Equations 

In the current and the following sections we will be interested in the equations for the 
bosonic components of the gravity multiplet, the graviton /i^,^ and the vector mode A^. 
We start with a ten-dimensional analysis of the linearized supergravity equations for the 
graviton excitations, valid for any solution on the baryonic branch, and proceed with a 
derivation of the equations for the vector field in the KS background in the section HI 

The traceless symmetric perturbation of the metric is described by the five-dimensional 
Klein-Gordon equation for a minimal scalar coupled to the background in the Einstein 
frame. A straightforward check [16j| shows that this property holds for the whole baryonic 
branch. 

Here we follow the notations of Papadopoulos and Tseytlin [8]. The functions A(t), 
p{t), x{t) and below are scalar functions from the PT ansatz depending on the radial 
variable t. In particular, A{t) is equivalent to the warp factor in the KS case e~'^^ = h^/"^. 
It should not be confused with the vector fiuctuation of the metric A = Ai dx*. In the 
Einstein frame the equation for the fiuctuation of the graviton 6 (ds^) = e~^^hfj_^dx'^dx'^ 
takes the form 

V + 2(i - $ + 2i) V - A^'e-^^-^^'-" V = 0. (11) 

where k'^ is the square of the 4-momentum and the over dots stand for t derivatives. The 
equation (fTTl) is precisely the Klein-Gordon equation for the minimal scalar in the the 
baryonic branch backgrounds including the KS point. 

To proceed to the explicit form of the equation (ITT!) for the KS background one chooses 



= h'/\ e^^+^x ^ 3 ^^^^^ ^ _ ^ ^^^^2^^ ^ 

e* = e*« = 1, e^^ = — (sinht cosht - t f^ h, 

16 



(12) 



where h{t) is the warp factor of the metric 



,,, , , (x cothx — l)(sinh 2a; — 2x)"'^''^ , , 

h{t)= dx^ (13) 

sinh X 



With these assignments the equation takes the familiar form [T7] 



•■ 8 sinh^t ■ 2 h{t) sinli^ t in, (-,a\ 

+ 3sinh2t-2t ~ (sinh2t-2t)2/3 = °- 
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In the last term we absorbed the numerical constants in the normalization of the momen- 
tum. It is also convenient to write the equation in the conventional Shroedinger form 



{-d^, + \/2(t)) V = 0, (15) 

with the effective potential V2(A;^,t) given by 

^ k'^h{t)smh^t 8 sinh^t 4 sinhtcosht 

^ ~ (sinh2t - 2t)2/3 ~ 9 (sinhtcosht - tf ^3 (sinhtcosht - t) ' 



4 Vector Mode 

To find a supergravity excitation that corresponds to the U{1)ti current JJ^, one should 
consider a special deformation along the angular direction d/dip of the T^'^ ~ S*^ x 5^, as 
was discussed in the section [2J We perturb the SU (2) x SU (2) invariant 1-form f llOl) 
in the following way: 

g^^g^ + 2/3(t)i, A = A^dx^ , (17) 

where A is a 1-form describing the vector mode and l3{t) is yet unknown function of t. 
Such a deformation leads to the following perturbation of the metric: 

ds^ ^ds^ + 2l{t)g^ ■ A, (18) 

where we introduced / = 2i3e~^^~^ for a latter convenience. This change of the metric 
will affect the Einstein equation as well as other equations of the type IIB supergravity. 
In particular one needs to modify the RR 5-form to preserve its self-duality: 

SF5 = -(3A A dg^ A d^^ + (3dA A g^ A dg^ + /Je^f+^Z^ *5 dA A dg^+ 

+ 2e-^%(3 - /3ir)e-=^P-"/2 *5 ^ A g\ (19) 

Here (3{t) is yet another function to be determined and K = AAe^^ in the PT notations 
[S]. This turns out to be a minimal ansatz required for the KS solution. One can show 
that there is no need to perturb the other type IIB fields if one is interested only in the 
four-dimensional transverse part of A. 

The ansatz so far contains the unknown functions /?, and 13 or /, which can be fixed by 
the equations of motion. The Bianchi identity provides us with the following equations: 

^^^g3p+x/2 ^ 2e-^%(3 - (3K)e-^P-^/^ *5 i = 0, (20) 

and a simple equation for the function /?, 

/9 = 0, or (3 = (3o. (21) 
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To find the function l3{t), or l{t), one should linearize the Einstein equation with the 
perturbation of the metric as in f|T8l) . The only nontrivial equation comes from the 6R^^ 
term. After certain simplifications one can write it in the form 

d^A^ + (2{i/l) + 6p + 3i + 2i) dtA^ - A;2e-2^-6p-i^+ 

+ ((VO + (V0(6p + 3i; + 2A) - 2A{Qp + x) - 2e-^2p-4xj _ 

= {Hi + Fi) - ^ e-^-^^ir + i e--ir^) i,. (22) 

In the KS background the square of the 3-forms is given by 

^ + 2 cosh^ t — sinh t cosh t + cosh^ t — 2 + cosh"^ t /„„x 

= H^ = ^e'^ 4 . (23) 

sinh t 



If one now writes the equation ( l20i) in components, taking into account (1211) and the 
transversality condition d^A^j_ = 0, 

d^A, + {6p + x + 2A)dtA^ - k''e~^^-^P-^A^+ 

+ {%~PA e-^2p-2-/;3p _ 2 e-i2p-4- j = q, (24) 

and compares it with the equation (l22l) . one will find that two equations coincide only for 

/3o = 1, and / = e"^. (25) 

Thus, the equation ^IM with the solution fl2^ describes the transverse vector excita- 
tion of the KS supergravity solution. For computation of the mass spectrum it is worth 
writing (!24|) in terms of the explicit solution (fT2|) . We obtain the equation 

dlA^ + V{t) dtA^ + Q(t) 4 = 0, (26) 

witlj^l 

^/ X 4 sinh^ t ^ ^ h 

Vit) = - - — - 2 cotht - - , (27) 

3 (sinh t cosh t — t) h 

k'^hsinh^t 8 sinh^t 2 hsink^t 

^' ~ ^ (sinh2t - 2tf^^ ^ 9 (sinh t cosh t - ty 3 (sinh t cosh t - t)/i ' ^ ' 

Again, one could write the above equation in the form f[T5l) with the new effective 
potential Vi(/c^, t), 

1- 1 . ^ k'^h sinhh ^ ^ ,2 1 (sinh 2t - 2t)^/3 

Vi = - V + - P^ - Q = ; TTTT^ -1 + 2 coth^ t + - ^ ^ — + 

2 4 (sinh 2t-2t) 2/3 4 /^sinh^t 

^ 3(sinh2t-2t)2/3(tcotht-l)2 ^ 2 tcotht-1 



4 /i2sinh^t 3 (sinh2t-2t)2/3/i 

2 (sinh 2t - 2t) (t coth t - I) coth t 
/i sinh^ t 



(29) 



Here we use the same momentum normalization as in the equation JT 
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Closing this section we notice that the equation fl26|l presented here coincides with the 
equation derived by Krasnitz in the UV hmit of the KS theory. The t ^ oo hmit of fl26l) 
is the same as the equation (4.30) of [1A\ with the assignment 

27 



and the change to the standard radial variable r = e 



i/3 



5 Super symmetry and 5d Approach 

In this section we compare our findings with the results obtained in the effective five- 
dimensional models of gauge/gravity correspondence [10] and show that the equations for 
the graviton and the vector mode are related by a Supersymmetric Quantum Mechanics 
transformation. This allows us to extend the equation for the vector mode to the baryonic 
branch. 

The authors of [TU] systematically study the 7^-symmetry invariant sector of fluctua- 
tions above the M = 2 backgrounds of the five-dimensional A/" = 8 gauged supergravity. 
Those also include the gravity multplet, i.e. the traceless four-dimensional metric fluctu- 
ation and the vector fluctuation, dual to the U{1)ti current. 

Although the KS solution truncated to five dimensions would correspond to a more 
general M = 2 supergravity theory , it is nevertheless interesting to compare the results 
of the two approaches. In fact, in both cases, the unbroken supersymmetry is A/" = 2 as 
we deal with the supergravity dual models of A/" = 1 gauge theories. Therefore the results 
based on the on-shell supersymmetry can be applicable in both cases. Indeed, we find 
that SQM transformations that relate the equations for the graviton and the vector mode 
in the case of the KS background coincide with the supergravity transformations used in 

In five-dimensional theories one can use the gauge freedom to recast the background 
metric into the kink form 

dsl = dq^ + e^^^^^rj^^dx^dx". (30) 

According to a general observation of [12j, the traceless graviton fluctuation /i^j, in 
five dimensions satisfies the equation for a scalar minimally coupled to the geometry (l30ll . 

(9j + 4r'9,-e-2^A;2) ^ q_ ^3^) 

Using the transformations of the effective M = 2 supergravity of [10] one can transform 
the graviton h^y into its superpartner - vector field 13^. As a result, the minimal scalar 
equation transforms into 

[d,e'^d,-k' + 2e'^^^^B, = Q. (32) 
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Here again /c is a 4-momentum. We are going to show that of fl26l) and i?^ are related 
by a simple field redefinition. 

The approach of ^lOj uses the superpotential, what can be problematic for the back- 
grounds from the baryonic branch (except for the KS solution) since the corresponding 
superpotentials are not known. Therefore there is a concern that the equations obtained 
for the KS may not be applicable for the outer branch. Nevertheless, we notice that the 
equation itself is ly- independent. This already suggests that it is actually valid for any 
background of the form (|30l) . Below we will give an argument based on supersymmetry 
that the equation ( |32l) can be applied to the whole baryonic branch. 

Let us first show that the equation (l32l) is the same as the equation (|26l) after an 
appropriate field redefinition. One can think of the metric fl30|) as an effective metric 
obtained by truncation of the ten dimensional theory with the metric (jl]) in the PT form, 
taken in the Einstein frame, 

dslo = (^e-^f-^dt^ + e^^r/^^dx^dx'^ + ^i^jdj/^dy^) 6"*/^ (33) 

The metric f l30p is then 

dsl = (e-s^-Mt^ + e^^r^^^dx^dx^) deti/3(^(5))g-4#/3 ^ 

(e-^P-Mt^ + e^^r^^^dx^dx^) e-2p+^e-^*/^ (34) 

what gives the following identification for the coordinate q and the function T{q): 

d ^ g4p+2#/3 d ^ 2T = 2A-2p + x--^. (35) 
dq at 3 

Hence the equations for the graviton in ten and five dimensions coincide, because they 
are just minimal scalar equations. 

The equation (!32|) in the PT notations takes the form 

d^B^ + (2p + X + 2i - ^ $) dtB, 
+ (^(4p+^<i>) (2i 

To compare this to fl26l) . derived in KS, set $ = 0. To match the kinetic terms in two 

equations one should redefine the field = e^^A^. After redefinition one gets 

d'^A^ + (6]3 + X + 2A) dtA^ - fc^e-^^-^^-^i^ + {2p (6i + 3x) + 2i + x) = 0, (37) 

which is precisely the equation fl26l) for the KS solution f|T2l) . 

We can further reduce the five- dimensional equations fl3Tl) and fl32l) to one dimension 
by taking the square of momentum k"^ to be the eigenvalue —im?. This will reduce the 
supersymmetry algebra to the Supersymmetric Quantum Mechanics with two differential 



2p + X - ^ $) + 2i - 2p + X - ^ $ ) 5^ = 0. (36) 
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operators Qi and Q2 that relate the solutions of the two equations (|3T1) and fl32|) . These 
operators realize the effective transformations of the supersymmetry algebra that was 
studied in \10\. Indeed, there are operators Qi and Q2, such that the equations 

QiQ2hf,u = -m^h^y and Q2QiB^ = -rn^B^ (38) 

coincide with the equation for the graviton (1141) and the equation for the vector mode 
f l26|) in the form (l36l) . It is easy to show that the operators that satisfy fl38|) are 



Qi = (a, + 2T') = e^^'+2*/3 (^^ + 2A-2p + x-^^ 



(39) 



and 

The operator Q2 is precisely the operator from (73) of [10] that realizes an A/" = 2 
supergravity transformation relating h^^ and B^. 

To get a more conventional representation of the SQM here, one can change the 
coordinates to dq = e~^du and bring the equations fl?I]) and fl5^ to the form f[T^ by 
redefining the wave functions and 5^. Let us define an operator 

Q-( (41) 



du + W 

with W = — 3T'/2, that acts on the vector made of redefined wave functions tph and ipB- 
According to the equations (13T|) and (l32l) the action of is as follows 

i'h \ _ 2 f i'h 



Therefore is analogous to the Hamiltonian of the SQM. Notice, however, that its 
eigenvalues are m^, not m, because Qi and Q2 correspond to the squares of the original 
supersymmetry transformations, i.e. Qi,Q2 are bosonic operators. 

We see now that the equation (l3T|l and (l32l) are related by supersymmetry transforma- 
tion for any background (!30|) . Since the minimal scalar equation describing the graviton is 
valid for the whole branch, the superpartner of the graviton (the transverse vector mode) 
satisfies the "superpartner" equation (152]) for any background from the baryonic branchjj 

We have calculated the spectrum of both equations numerically for the backgrounds 
along the baryonic branch. Since the equations for the superpartners are significantly 
different the discrepancy between the masses can be used as an error estimate of the 
numerical method used in the calculation. 



In general, there is a family of equations like l|32]) that are related to (|3T]l by a supersymmetry 
transformation. Indeed, for a given W from (|¥T|) . any W that satisfies + W' = + W gives rise 



to such an equation through l|42)) . Nevertheless, the equation (|32)) is uniquely specified by a requirement 
that the effective potential Vi is singular at i = 0. This is true because Vi is singular in the KS case (|29)) 
and hence should be singular everywhere on the branch by continuity. 
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6 Numerical Analysis 



In this section we present the results of the numerical studies of bound state spectra 
for the baryonic branch backgrounds. In our computations we will rely on the shooting 
technique. The spectrum of the minimal scalar equation f|T^ in the KS background was 
also studied numerically in [161 [IS [IS] while the analytical approximation was employed 
in [13]. 

We start by comparing the KS spectra of the equations for graviton f[T^ and vector 
mode fl2Bl) . Two fluctuations are related by supersymmetry and thus their masses should 
be the same. The spectrum is presented in the table [1] The eigenvalues match with those 
obtained by Krasnitz [17] with the WKB approximation. Comparing the numeric values 
of the masses of the spin-2 and vector particles in the table [1] one could estimate the error 
of the shooting technique in the KS case to be around 0.1%. 



n 


1 


2 


3 


4 


5 


6 


7 


8 


9 


Graviton 


1.764 


4.002 


7.143 


11.19 


16.16 


22.03 


28.83 


36.54 


45.16 


Vector Mode 


1.762 


3.999 


7.136 


11.18 


16.12 


22.01 


28.80 


36.50 


45.12 



Table 1: The spectrum of for the gravity multiplet 



First few (up to ten) values of in the KS spectrum can be approximated with a 
good accuracy by a quadratic fit 

= 0.46^2 + 0.86^ + 0.46, n= 1,2,3,... (43) 

We present the results of the fit and the masses on the figure [T](a). It is interesting that the 
fit 0431) is close to the spectrum even for small n. The fitting formula is proportional 
to (n + no)2, where Uq is close to one. This is consistent with the approximation of [13] , 
where the eigenvalues were matched to zeroes of the Bessel functions, ubiquitous in the 
conical geometry. A similar result was obtained in [20] for the GPPZ [21] flow, where the 
exact spectrum was proportional to (n + 1)^. 

The fit fH5]l was found by minimizing the sum 

N 

\m\ - [c^r? + cin + cq) f (44) 

n=l 

for the few first states = 5, . . . , 10. With more points taken into account the least 
square fit would increase the accuracy of the highest coefficient by the price of a larger 
deviation from for small n. We found C2 to be ~ 0.459 in the KS case. This number 
is in good agreement with the universal coefficient obtained by Berg, Haack and Miick in 
[T9] . In their normalization the coefficient takes value (3/4)^/^/1(0) C2 — 0.27. 

Remarkably, the coefficient C2 does not depend on the details of the effective potential, 
but rather encodes information about the background geometry, namely, the combination 
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(a) 



(b) 



Figure 1: (a) Values of for the graviton multiplet in KS for different quantum numbers n. (b) Ex- 
tension of the spectrum to the baryonic branch parameterized by U . 



g^^Qtti which arises from the Laplace operator in five dimensions. Indeed, the WKB 



approach, apphed in [TT], gives 

/ dt ^-V2{t) 

Jo 



TT + (ra — l)7r. 



(45) 



where V2{t*) = 0. In the KS case V2 is given by (fT6|) . For large n, and consequently large 
TUn, the ^^-independent term in V2 can be dropped and we obtain an analytical expression 
for C2 in the KS case 

\/h sinh t 



C2 = TT 



dt 



(sinh 2t-2t)V3 



0.460 . 



(46) 



Let us choose the coordinate U, introduced in [5] , to parameterize the baryonic branch. 
To estimate the scale of the spectrum for a non-KS background we rewrite the potential 
f|T6|) in terms of the PT ansatz [8J, substituting k"^ for its eigenvalue — m^: 

[a cosh t + lf + 2a^ sinh^ t 



V2{m\t) 



+ 



2a cosh t 



-39 



where a{t) is another function from the PT ansatz [8|, e^^ 



(47) 
2a cosh t, and 



Although we cannot find the spectrum of analytically, we can estimate how it 
scales with the parameter U when we are significantly far from the origin of the branch. 
We start our analysis with the m^-independent part of V2, which only slightly varies as 
we increase U. Indeed, its leading UV (t 00) asymptotic is [/-independent: 



^2(0, t) = - 



4 (5 - 2t) 



6 



-At 12, 



+ 



(4J 



12 




-300 -200 -100 100 200 300 



U 

(a) 

Figure 2: (a) C2-coefficient as the function of 
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[/-parameter, (b) logC2 as the function of log f/. 



and ¥2(0,1) varies within a small range in the IR (t = 0): 

V2{0,t) = ^-^ai~0 + 0{t'). (49) 

Here we remind that ^{U) G (1/6. ..5/6) [4j is a function of U, which can also be used 
to parameterize the branch. It varies within the specified limits, and the point = 1/2 
corresponds to the KS solution. Hence 1^2 (0, 0) = 2/5 for KS and ¥2(0, 0) approaches 1/3 
for large U. The ¥2(0, t) is monotonic and therefore it can be approximated by a constant 
in the analysis below. 

Unlike V2{0,t), the mass-dependent component 171"^ e~'^^~^^v~^ significantly depends on 
U. It monotonically changes from a finite value at zero to the zero value at infinitjl^ 



-2A+a;„ 



.1 21/33 



e -'-V ^ = (4t_i)e-2*/3 + ___ (50) 
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In general, the value at zero is a complicated function of U, ^{U) and = ^{U,t = 0). 
It can be simplified in the large U range by substituting the limiting value = 5/6 and 
expressing $0 in terms of U and ^ [5]: e*" ~ 2^^'^3~^^^U~'^^^. This gives 



.1 21/33 



2U 



l_e2*o(i + ^ + ^ + 
9 135 



(51) 



The normalized solution to the equation (1151) exist only if V2 < at the origin, which 
suggests that scales at least as U for large U. We can try to be more precise using 
the semiclassical approximation and express the n-th mass through the integral over the 



^Here we use the normahzation of the warped factor introduced by Krasntitz [17j . what resuhs in 

-2A ^ 21/33^(6-2* „ 
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V2 < region, as we did above in fH5l) . This integral can be roughly approximated 
as a/— V2(0)t* ~ mt*f/~^/^. The main complication is to estimate t*. Since e^*° from 
f lSTj) is small, the perturbative expansion flSTl) suggests that t* increases with U until a 
point, where fl^T]) is no longer reliable. At the same time the large t asymptotic fISUl) 
is [/-independent, what suggests that for large U the value of t* approaches a constant. 
Therefore we expect ~ [/n^ for sufficiently large U. 

Numerical studies of the graviton multiplet spectrum on the baryonic branch shows 
the pattern depicted in the figure [TJ^b). Calculations confirm that the leading coefficient C2 
grows as f/", where a approaches 1 for large U (figure [2]). As a final touch, we collect in the 
table [2] the known evidence about the U scaling parameter a for some non-perturbative 
objects on the baryonic branch. 



SUSY D5 Baryonic Condensate Fundamental String 


Glueballs 


D3,D3 


a a <1 1/4 


1/2 


5/4 



Table 2: Scale behavior for large U: T ^ U' 



7 Discussion 

In this work we present the equations describing the bosonic degrees of freedom of the 
gravity multiplet for the KS and the baryonic branch backgrounds. The equations were 
derived by a linearization of the ten-dimensional type IIB supergravity equations. The 
traceless graviton from the gravity multiplet satisfies the equation for a scalar minimally 
coupled to the background f|TT]) . The vector mode of the gravity multiplet dual to the 
U{l)'ji current satisfies the equation fl2^ in the KS background. Its generalization to 
the baryonic branch fl5B]) is found by matching it to the equation fl52]) derived in the 
five-dimensional approach in [10]. This result is supported by the supersymmetry trans- 
formation that relates the wave functions of these fiuctuations. 

The mass spectrum of the gravity multiplet for the KS background can be found in the 
table dJ This spectrum can be approximated with a good accuracy by a simple quadratic 
formula fH5]) . which is approximately 

~ 0.46(n + 1)2. (52) 

This simple complete square form does not hold along the baryonic branch, although the 
spectrum can be well approximated by a general quadratic formula C2n^ + Cin + cq. 

In this work we did not study the spectrum of the anomaly mutliplet S* of ([1]), which 
contains the fiuctuation of the trace of the metric hj^ and the longitudinal part of the 
vector fiuctuation a ([5]). The main complication is that these fiuctuations do not decouple 
from the other supergravity fields. Recently the fiuctuation of the metric hj^ was consid- 
ered as a part of 7-particle system by Berg, Haack and Miick in [13 [H] . They found the 
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resulting spectrum of the system, but the individual mass towers were not identified with 
the glueballs. 

Based on the similarities between the spectrum of the gravity multiplet in the KS fl52l) 
and GPPZ backgrounds, where = 4L~^(?t, + 1)^, one can assume that some features of 
the spectra for certain glueballs do not crucially depend on the details of the background. 
Based on the exact result of the GPPZ case calculation for the mass spectrum of the 
anomaly multiplet S, ml = AL-^{n + l)(n + 2) pB[20l[22] , one can guess the answer for 
the KS case. In the units of BHM the approximate formula reads 

~ 0.27(n+ l)(r2 + 2), n = 1,2,3,... 

This is in fact close to the lightest of the seven towers of BHM, given by the empirical 
formula 

0.271^2 + 0.774n + 0.562. 

It would be interesting to confirm the matching between the trace of the metric and the 
lowest tower of the 7-particle system with a more rigorous approach. 

We find it intriguing that the spectrum of [19] contains only two states that look 
degenerate. Given that we are dealing with massive states of the = 1 system we would 
expect all the states to be degenerate. This might signify that the numerical method used 
in [19] alter the mass degeneracy because of a numerical error. Another possibility is that 
the superpartners of the glueballs in question are not the part of the 7-particle system 
and can not be captured by the fluctuations of the PT ansatz. The clarification of the 
magnitude of the numerical error is also important to check another finding of [IH] - the 
significant deviation of the spectrum from the quadratic behavior for few lowest values of 
the quantum number n. 

We are grateful to M. Berg, O. DeWolfe, M. R. Douglas, D. Z. Freedman, I. R. Kle- 
banov, A. Konechny, G. Moore and N. Seiberg for very useful comments and stimulating 
discussions. The work was partly supported by Russian Federal Agency of Atomic Energy; 
by the Council of the President of the Russian Federation for Support of Leading Scien- 
tific Schools NSh-8004.2006.2, and by the grants NSF PHY-0243680, RFBR 07-02-00878 
(A.D.) and DOE DE-FG02-96ER40949, RFBR 07-02-01161 (D.M.). 
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